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Introduction 


I-'iiut(‘-('lriu('nt  uiudcls  arc  used  ('xl('iisivelv  l)v  the  rotorcraft  iiidiistrv  for  deteriiiiiiiiiy,  the  vil)ratorv 
characteristics  of  ia'licopK'r  ai'draines.  These  airfratii<‘  itiodels  are  r\’picall\'  large.  ii.suall\'  coiitaiiiiiig  .se\-eral 
thousand  tniite  eltuiKuUs.  and  their  definition  involves  eonsidertible  effort.  Hesnlts  olttained  by  using  t  lu'se 
models  tisuallv  do  not  correhite  well  with  ex])ernnent al  d;ita.  particularly  at  the  higher  freciuencies  that  are 
near  the  predominant  rotor  excittuion  f'r('(|nencies  (refs.  1  through  4).  .Although  there  art'  severtil  possible 
retison.s  foi'  the  lack  of  correhit  ion.  oik'  that  has  rt'cently  attracted  attention  is  the  accurac\'  of  the  elements 

t  henist'lvcs. 

-Most  tinit t'-elemeiit  codes  in  use  lodtiv  ('nii>loy  wliat  ea'C  eonniKnily  referred  to  as  h-version  elements.  .An 
hwersion  I'leinent  list's  iixeil-order  shape  functions  to  n'late  the  discrett'  nodal  displact'inent  s  of  t  lit'  element 
to  the  continiiou.s  displacements  within  the  ek'inent.  .Mesh  size  normally  controls  the  accuracy  of  the  rt'sults 
!ibt. lined  Using  li-version  (  if'ments.  .As  the  niiinbci  ol  elements  is  increased,  the  mesh  size  decretises  and  the 
result.'  it  i.d  it/  .  , ig;'.  ic.s;.  .'or  coiu’t'rgi'nce  usually  involve  solving  a  series  of  jiroblt-ms  with  snccessivt'lc- 
refined  meshes  until  convergence  is  indicated.  .Such  convergence  checks  are  often  used  in  studii's  which  emplo\‘ 
simdl  mode!-'  To  che'  k  the  linio'-i'lement  model  of  a  hirge  structure  (such  as  an  airframe)  for  convergence  b\' 
mesh  relinement  would  be  ;i  formidable  ttisk  not  likely  to  be  done  in  practice.  So-called  p-version  I'li'ments  th.tt 
hiivi'  recentlv  gained  titti'ution  ha\'e  the  [lotential  for  both  improving  accuracy  and  simplifying  convergence 
checks. 

The  p-v('r.sion  cl(>in('nt  is  differt'ut  from  the  h-version  el('ment  in  that  tht'  ust'r  may  si'li'ct  tht'  ordt'r  of  tht' 
shapt'  function  defining  the  disjiiact'inent  behavior  of  tht'  elemt'tit.  In  gt'iiertd.  the  higher  tht'  ortli'r  ttf  the  shttpe 
function,  the  mort'  ticcurati'  the  rt'sults.  Increasitig  tht'  tirtlt'r  t)f  the  shap<'  functions  is  analogtms  to  incretising 
the  number  of  t'lt'mciits  for  ;i  com[)tirabli'  h-vt'i'sitni  intitlt'l.  This  featnrt'  becomt's  a  distinct  tiilvantage  in 
convergt'iict'  checks  of  large  fini'e-t'lement  modt'ls  bt'canst'  raising  tht'  ttrili'r  of  thi'  sh;»i)i'  functions  would  be 
much  I'asit'f  tlnm  increasing  the  numlit'r  of  ek'int'nts  in  tht'  thita  input  fik'. 

Tilt'  ativanfage  of  using  [i-vt'rsion  ciemetits  is  lost  if  the  element  tiot'.s  jiof  clost'ly  upproxiiiintr  flit'  gt'otnefrv 
of  thf'  structure  being  mtnlek'd.  For  t'xtunple.  a  nniftirm  [nvt'rsitm  bttam  element  vvoukl  ptist'  nt)  athantage 
over  ;i  tmifortn  h-vt'rsion  t'k'inent  in  t'ither  accuracy  tir  t'ast'  t)f  u.se  in  modeling  a  tapt'it'il  bt'tun.  Howevt'r. 
a  p-version  beam  I'lement  whose  cross-st'ctitmal  tiimt'iisittns  vary  linetu'iy  with  k'ligth  would  allow  ftir  a  mtire 
accuratt'  represt'utation  of  a  tapereil  bt'am  titan  ettuki  bt'  achit'xetl  with  a  uniform  bt'am  t'k'ini'nt.  Bt'cause 
many  substructurt's  in  a  typical  airframe  ctuisist  ttf  beams  w'hose  cr-tss  sections  vary  aktng  tht'ir  k'ligth.  tht'  use 
of  p-vt'fsion  elements  offers  an  ativanfage  ttver  li-versittn  elements  ttniy  if  the  [t-versittn  elemi'iits  can  accuratt'ly 
motlel  tht'  physical  structure  without  resorting  tt)  exet'ssive  mesh  rt'fini-mt'iit . 

The  Langk'V  Ifesi'arch  Center  has  nntk'rway  a  rt'search  structural  tiynamics  jtrogram  which  is  invi'st igat ing 
finite-ek'iiK'nt  modt'ling  techiTuiues  fttr  lu'licttpter  vibratittns  analystts.  Tnis  activity  has  itk'ntiiii'tl  tht'  ni't'd 
for  a  small  rt'search  code  containing  a  limitetl  library  ttf  p-vt'rsittn  bt'am  anil  slit'll  ek'int'nts.  'Ttt  inet't  tht' 
ri'tjuirements  of  a  research  tool,  these  elements  must  be  well  dttcmnt'ntetl.  t'asily  mtttlifii'tl.  anil  tailoretl  to 
the  net'ils  of  airframe  vibratittns  analyses.  \\’ttrk  lias  Itet'ii  initiatetl  at  Langk'y  tti  devt'kip  such  a  cotie.  Inititd 
attt'ution  is  bi'ing  given  to  the  beam  t'lement.  ,A  survey  ttf  the  finite-element  litt'ratiiri'  rt'vt'aletl  ntt  cttnnni'rcially 
available  codt's  with  it-vt'rsion  beam  elemt'iits  anti  tinly  a  few  research  ctitit's  having  such  elements,  .\tine  of  the 
latter  coik's  appear  to  bt'  suitable  for  airfratne  vibratittns  analyses  wttrk.  Httwt'vi'r.  tht'rt'  has  bt'en  significant 
work  done  in  tht-  ari'a  of  taitt-retl  h-versittn  beam  elements  w'hich  is  rek'vant  to  tht'  bt'am  t'lement  that  is  thi' 
subject  of  this  jiajier.  Pt'rtinent  work  in  this  area  is  summarizt'tl. 

Tapered  beam  finite  ek'tnents  basetl  on  cubic  pttlynttmial  shajtt'  functions  havt'  bt't'ii  iintler  invi'st  igat  ion 
for  many  years  (refs.  5  anil  fi).  Highi'r  oriler,  tajtereil  beam  eli'ini'iits  havt'  also  bet'ii  devt'kipi'il.  'Thomas  anil 
Dokiimaci  (ref.  7)  preseiitt'il  two  ta[ti'reil  beam  elt'U..  .ps  usit'g  t;uinti.  itol v nomials.  One  bi'atn  I'k'inetit  had 
*"•<'  i'lti  .iiai  iM'O'  v\  nci'i  as  tile  oilii'f  bt'am  element  tlefinetl  tht'  st'conil  derivative  at  t'ach  I'lid  of  thi'  bt'am  as 
a  degree  of  fft'cilom.  Inferelement  continuity  was  then  enforcetl  for  ilisplacement  and  thi'  first  two  dt'rivat ivt's. 
which  fornieil  a  C'^-type  ek'tnent.  .A  tapt'rt'il  beam  i'lement  using  st'vent h-tk'gri'i'  polynomial  shapt'  functions 
was  derived  by  To  in  reference  IS.  Kach  noik'  hail  four  ilegrt'i's  of  fret'tlom  (ilisplaci'ini'iit  ami  tht'  first  three 
derivatives)  on  which  intert'li'tni'iit  continuity  was  impost'll.  It  should  bt'  poititetl  out  that  these  higher  order 


are  ’int  p-vrsiou  ('leniciits  Ix’causc  the  onicT  of  tlio  shape  functions  was  Hx('<l  when  th('  element  was 
derived.  In  another  paper.  To  (ref.  9)  d('veloi)(‘d  a  tapered  beam  element  incorporating  shear  deformation  and 
rntarv  inertia  but  ba.sed  tlu'  element  formulation  on  cubic  shape  functions.  .Mon'  recently,  a  tapered  p-version 
beam  el.  lueiir  was  dt'velopi'd  by  Hodges,  Hopkins.  Knnz.  and  Hinnant  (ref.  10)  s[)ecifically  for  modeling  rotor 
blades.  To  correctlv  niodt'l  a  spinning  rotor  blade',  the  e'lenu'iit  inclneles  the'  nonline'ar  e'ffe'cts  of  large  nodal 
displace'iiu'iits  and  rotation>,  at'rodynamics.  and  inertial  rotation.  Such  a  complex  element  is  not  needed  for 
modeling  airframe'  st nu  t  lire's.  .\s  a  practie'al  matter,  elements  u.sed  feir  the  analysis  e)f  airframe's  sheeulel  be' 
relati\i'l\'  simple'  be'e  ause'  of  the-  large'  numbe'r  e)f  elements  normally  re'einired  te)  meeelel  the  strue-tnre  aele'einate'ly. 

rill'  objee  tive'  of  this  papi'r  is  to  presi'iit  the  ele'rivation  of  a  tapere'd,  i)-ve'rsion  be'am  e'le'ine'iit  for  use  in 
dviiamii-  analysi's  of  ge'iu'ral  strue-tural  systems.  This  element  leature'S  hie'rare-hie-al  shape-  fune-tieens  whiedi  allow 
higher  orde-r  analyse-s  to  use'  ele'ine'iit  matrie'i's  e-stablishe-el  for  lower  oriler  analysi's.  .Appropriate-  orthogonal 
re'lations  for  the'  shape-  fimitions  are-  i-mploye'il  to  avoiel  ill-e-onelitione-el  matrices  anel  re-dui'e  the-  number  of 
non/ero  te-rms.  Eli-mi'iit  matriee-s  are-  i-xplie-itly  forme-el.  thus  eliminating  the  ne-eel  for  numerieal  eiuaelrature 
ri'siiltiiig  in  a  simple-r  impli'ini'iitat ion  anel  a  ri'dne-tion  in  ronneloff  e-rror. 
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I'oe'flie  ients  of  i)olynomial  re-pre-se-nt ing  cross-sectional  property 
modulus  of  eiasticity 

etli  axial,  laii'ial  be'iieiing.  and  torsional  freeiue-ncy. 
ri'spe'ctively,  Hz 

slie-ar  moelulus 

number  of  internaf  degrees  of  freedom 

are-a  mome'nt  of  inertia  about  .V-.  V-.  anel  Z-axis.  respectively 

torsional  stiffne-.ss  e'onstant 

stiffness  matrix 

e-xte-rnal  stiffness  snbmatrix 

e-xte-rnal-internal  stiffne'ss  e'onpling  submatriee-  assoe’iateei  with  ». 
c.  It',  anel  reespee-tively 

internal  stiffni-ss  submatriee  assoe  iati-ii  with  it.  c.  tr.  anel  (K 
re-spee'tively 

li'iigth  of  beam 

mass  matrix 

exte-rnal  ma.ss  snbmatrix 

I'Xte-rnal-inte-rnal  ma,ss  eonpling  submatriee- a.ssoe'iati-el  with  n.  r, 
ir.  anel  f).  re-spi'e  tive'ly 

inte-rnal  mass  snl)matrie'e  a.ssociate'ei  with  ii.  r.  ir.  anel  fJ. 
ri-spe-ctively 

total  numbi-r  of  lie-gri'es  of  fre'e-elom 
numbe-r  of  i-lemi-nts 

dh  shape-  function  associate-el  with  it.  r.  if.  anel  0.  ri-spi-et ively 
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all  shape  function  of  typo 
/th  shape  function  of  type 

nuniiier  of  shape  function  associated  with  u.  v.  w,  and  0. 
respectively 

qtiadratic  polynomial  representing  generic  cross-sectiona' 
property 

ciuartic  polynomial  representing  generic  cross-scctional  property 

dh  discrete  degree  of  freedom  associated  with  n.  v.  tv.  and  0. 
respectively 

matrix  function  defined  by  /q  P4(x)NfNj  dr 
matrix  function  defined  by  Jq  P4(x)( dx 
matrix  function  defined  by  Jq  P2{x)N^ Nj  dx 
matrix  function  defined  by  /q  P4(x)(Af^ dx 
matrix  function  defined  by  /q  P4{x){N^  )"(Xj)"  dx 
kinetic  energy 
time,  sec 
strain  energy 

displacement  along  A'-.  Y-.  and  Z-axis.  respectively 

rectangular  axis  system  with  origin  at  one  end  of  beam: 

AT-axis  along  neutral  axis  of  beam,  and  V'-  and  Z-axes  oriented 
parallel  to  principal  axes  of  beam 

independent  variable  along  A’-.  Y-.  and  Z-axis.  respectively 
coefficients  in  generating  equation  for  jV-' 
coefficients  in  generating  equation  for 
first  variation  of  (  ) 

r  0  j) 

ll  {i=j) 

j  -  I  (f  =  2  or  4) 

1  (Otherwise) 

rotation  about  A'-axis  (i.e..  tonsion) 
ma.ss  density 

fth  axial,  lateral  bending,  and  torsional  frequency.  res{)ectively. 
rad/.sec 


f’rimes  to  a  svmbol  denote  a  derivative  with  respect  to  x.  A  dot  over  a  symbol  denotes  a  derivative  with 
resp<'ct  to  time. 
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Basic  Mathematical  Formulation 

riu'  hcaiu  is  dcrua'ii  assuiniiif^  that  tlio  hcam  hchavcs  ki''''inatically  like  a  Bernoulli- Eiuk'r 

heaiii.  Horar\'  inertia  ('Meets  ar('  ineluded:  li()\v<'V('r.  shear  Hexihility  is  not  inehuh'd  in  the  prf'sent  formulation. 
1  he  eros.'-seer ional  aia'a  ran  \-ary  as  a  ([uadratie  polynomial  along  the  k'ugth  of  the  heam.  and  the  aia'a  moment 
ol  inertia  can  \;u\'  a,--  a  (luanie  ixilynomial.  Fhe  ('h'luent  ;ixis  syst('m  \.\.Z  is  oriented  .so  that  the  A’-axis  is 
along  the  neutrtd  a'ds  of  thf'  midefornu'd  Ix'iim,  and  tht'  )  -  and  Z-axes  art'  parallel  to  tlu'  principal  axes  of  tlie 
eros>  section.  '  .Sc('  Hg.  1.)  This  orientation  rt'stihs  in  a  zero  cross  jtroduct  of  inertia  =  0).  The  continuous 

( iisplacements  ii.  c,  u\  and  d  tire  tissumt'd  to  Ix'  oi.Iy  a  function  of./-  and  time. 

Hiuaiiton's  principh'  is  us('d  in  tlx'  deriviition  of  this  heam  element.  .Assuming  only  consf'rvat ivt'  forces. 
Hamilton's  principh'  is  statt'd  as  (ref.  11) 


t , 


f'iT  —  ri  lit  =  (t 


(1) 


.It 


wlx're  and  hj  ri'pia'Sf'ut  arhitrar\'  tinu's  at  which  the  state*  of  thi'  systi'm  is  known,  and  T  and  T  are  the 
kiiK'tic  and  strain  ('lu'rgif's.  ix'spe'ct i\-('ly.  The  kiiu'tic  energy  of  tlx*  heam  is  found  in  la'h're'net'  12  as 
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[  I'l  r 

--  /  /)  An~  -h  A>''  +  Air~  +  I  \  -I-  I  /  y(i'' I)')' ( u’' )' 

2  ./(I  L 


(l.r 


(2) 


and  the  strtiin  ('iK'ruv  follows  sin.iila.rl'.'  frotn  r<‘ferenc<'  13  ;is 


I'  = 


EAiu')-  +  Elyzlr")-  +  Elyyiw'']^  +  d.r 


"A 


j/,2 


(3) 


The  continuous  problem  associatt'd  with  the'  continuous  displact'nH'iits  u.  c.  le.  and  0  is  discretized  by 
introducing  discrete  de'gre'es  of  freedom  tj,  which  art'  related  to  the  continuous  displact'iiK'Uts  according  to 


II 

(■1) 

/>, 

,.  =  X^.V 

1=1 

Id', 

(a) 

II 

N/,"' 

(0) 

II 

f  h 

(7) 

.Substituting  ('([nations  (1)  through  (7)  into  e(|natiotis  (2)  and  (3)  expresses  the  kinetic  and  strain  ('ix'rgit's 
in  t('rms  of  discrete'  degrees  of  fre-edom.  The  a[)plic;iti()n  of  Hamilton's  jtrinciple  (('([.  (1))  the'ii  leads  to  the 
ieh'iit ificat ion  of  the  following  nonzero  terms  which  a[)pear  in  the  eleme'rit  mass  and  stiffness  maiiTe's: 


■'C  = 


/' 

./() 


'.v; 


d.r 


m:[j  =  £p[..tv;'.v;' +  A7'i'(.v;')' 

m;[j  =  |^([T.\7'■.v;'  +  />■p(.v;'•)'(.v;'j']  d.. 


(S) 

('.)) 

(10) 
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Mf  ,  =  pIxxNfNf  ds 

(11) 

'*-i 

7'"-' 

II 

(1'2) 

K[j  =  Eizz(^fy'(^jy'  dx 

(13) 

Kfj  =  l‘  EIyy(Nff(Nfy'  dx 

(14) 

Kfj  =  G.J(Nfy(Nfy  dx 

(15) 

Shape  Functions 

The  discrete  degrees  of  freedom  are  divided  into  two  sets,  external  and  internal.  The  12  external  degrees 
of  freedom,  which  are  depicted  in  figure  2.  correspond  to  the  usual  definition  of  the  physical  nodal  degrees 
of  freedom  for  a  beam  finite  element  (ref.  14).  The  internal  degrees  of  freedom  have  no  physical  significance' 
but  are  simply  the  coefficients  of  the  higher  order  shape  functions.  These  internal  degrees  of  freedom  serve' 
te)  ine'rease  the  ae’ceiracy  of  the  transformation  from  the  discrete  problem  having  a  finite'  number  of  elegrees  eef 
freedom  to  the  e'ontinnous  problem  having  an  infinite  number  of  degrees  e)f  free'elom.  The  number  of  internal 
elegrees  e)f  freeelom  used  in  the  beam  element  is  specified  by  the  data  input  file  and  may  vary  from  zero  to 
theoretically  infinite.  Setting  P„  =  Pg  =  2  and  P,r  =  P,,.  =  4  in  eepiations  (4)  through  (7)  will  lead  to  the 
classical  12-degree-of-freedom  beam  element.  Incre'asing  any  of  these  P's  adds  internal  degrees  of  freedom  to 
the  element.  The  shape  functions  Nj  are  usually  taken  to  be  polynomials  although  in  theory  they  can  be  any 
set  of  functions. 

Shape  Functions  for  u  and  6 

The  shape  functions  and  Nf  are  identical  and  have  C**-type  continuity.  That  is.  continuity  is  enforced 
across  element  boundaries,  but  continuity  of  the  derivatives  across  element  boundaries  is  not  enforced.  Shapt' 
functions  satisfying  continuity  will  be  denoted  by  Nf  (i.e..  Nf  =  Nf  =  Nf).  The  first  two  shape  functions 
in  this  set  are  the  well-known  linear  imlynomials  (ref.  15) 


A’l”  = 


-I-  1 


(lb) 


Nf  =  j  (17) 

The  higher  order  C'^-type  shape  functions  used  herein  were  deriverl  subject  to  two  reciuirements;  First,  the 
continuity  is  enforced  by  restricting  the  higher  order  shape  functions  to  be  zero  at  the  element  boundaries, 
and  second,  the  set  of  higher  order  polynomials  must  be  orthogonal  with  respect  to  their  first  derivative. 
Orthogonality  of  the  first  derivative  was  chosen  over  polynomial  orthogonality  Irecause  the  ('lenient  mass  and 
stiffness  matrices  obtained  by  reepnring  first-derivative  orthogonality  contain  fewer  nonzero  terms:  thus  explicit 
integration  is  facilitated.  This  also  results  in  matrices  which  are  better  conditioiu'd  than  those  obtaim'd  from 
orthogonal  shape  functions.  These  requirements  are  expre.ssed  mathematically  by  the  following  ecpiations: 


Nf{0)  =  0 
Nf{l)  =  0 


(i  >  ;i) 

(i  >  3) 

('>3../  >3) 


(18) 

(19) 

(20) 
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riu'  first  three  higher  ord^'r  C’^’-type  shaite  fiiiietioiis  in  the  sot  dofinod  by  equations  (18)  tlirougli  (20)  are 


•vJ  =  ^  f )  (-) 

In  general,  once  .V^’  is  known,  the  i  +  1  shape  function  can  be  found  from  the  recursive  forinida 

=  (of  +  (i-  >  3)  (21) 

where 


(21  -  3)v/2/  -  1 

{i  >  3) 

(25) 

i\'2i  —  3 

.If  = 

(41  -6)v/2r=T 
i/v/21  -3 

(i>3) 

(20) 

(1- 3)^27^ 

{i  >  3) 

(27) 

lV'21  -  5 

Shape  Functions  for  c  and  w 

The  shape  functions  X-'  and  iV."'  require  C^-type  continuity,  meaning  tliat  both  the  functions  and  their 
first  derivatives  must  be  continuous  across  element  boundaries.  The  shape  functions  N-'  and  A^“'  are  identical 
except  that  A'.]''  =  find  Nf  —  -N^  to  ensure  that  the  discrete  rotational  degrees  of  freedom  %  and 
have  the  sense  indicated  in  fi„ure  2.  Shape  functions  satisfying  continuity  will  be  denoted  by  (i.e.. 
.\f  —  A(()A'“'  =  A',^).  The  first  four  C^-type  functions  for  the  beam  element  are  (ref.  15) 


^3  ,w,2 

jsj}  =  2- _ 3—  4-  1 

.Vi  /  ^3  3  ^2  +  ^ 


X 


N^-^-2-+x 


r3  t2 

yvl  =  —  -  — 

^  /2  I 


(281 

(29) 

(30) 

(31) 


The  derivation  of  the  higher  order  C^-type  shape  functions  employed  in  this  paper  is  based  on  the  same 
|)hilosophy  as  the  derivation  of  the  C^-type  shape  functions.  To  ensure  continuity  the  higher  order  ,shai)e 
functions  must  have  zero  slope  and  displacement  at  the  element  boundaries.  Consistent  with  the  i)r('vious 
discussion  on  the  orthogonality  properties  of  the  C^-type  shape  functions,  the  C’-type  shape  functions  are 
rcfjuired  to  be  orthogonal  in  their  second  derivative.  The  specific  requirements  are 

Af/(0)  =  0  (i  >  5)  (32) 
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(i  >  5) 

(i  >  5) 
('  >  5) 
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(43) 


Matrix  Functions 

Inspection  of  equations  (2)  ami  (3)  indicates  that  there  are  eight  cross-sectional  properties  associated  with 
a  beam  element.  Two  of  them.  pA  and  EA.  are  represented  by  the  quadratic  polynomial 

P2(x)  =  02^  +  niJT  +  qq  (44) 

The  remaining  six  properties,  pl\x-  pWy ■  P^ZZ-  GJ.  Elyy.  and  EI^z-  ^^re  represented  by  the  <}uartic 
polynomial 

<  *1  o 

Pi(x)  =  a^x  +  a'iX'  +  (i2X  +  <i\x  +  no  (45) 

Because  these  polynomials  can  represent  any  of  the  aforementioned  cross-sectional  properties  by  substituting 
in  the  appropriate  values  for  the  coefficients  no  ecjuations  (44)  and  (45)  are  referred  to  herein  as  "generic 
cross-sectional  propjerty  p)olynomials."  Note  that  P2(x)  is  a  subset  of  P.i{x).  and,  therefore,  the  coefficients  uq 
through  (12  refer  to  both  polynomials.  The  context  in  which  the  coefficients  are  used  will  clearly  sj)ecify  which 
generic  cross-sectional  property  p)olynomial  is  being  referenced. 

Substituting  the  generic  cros.s-sectional  property  polynomials  for  the  actual  cros,s-sectionaI  prop)erties  in 
equations  (8)  through  (15)  and  substituting  generic  or  .V'  .shap)e  functions  for  the  displacement-sp)ecific 
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slia])!'  functiotis.  flie  10  ti'rins  in  ('([nations  (?<)  thronsh  (15)  arc  rt'dnccd  to  5  nni(|n('  terms.  Tlu'sc  five  nni(|ne 
terms  aia'  referred  to  liert'iti  as  "matrix  fnnet ions.”  and  have  the  follovvinf>,  definitions: 


Sl\j  =  l‘  ,U- 

dr 

S][j  =  l‘  r2(r)\Y.\;  dr 
•S’/j’  =  l‘ <lr 
=  I'  Pi(r)[.\Y)"(.\^f  dr 


(40) 

(47) 

(48) 


(49) 


(50) 


The  matrix  fnnetions  are  functions  of  a  cross-sectional  [trojx'rty  which  is  shown  as  an  argument  in 
[tarent iies('s  when  a[)|)r()[)ri;ite.  for  example.  5j’^(p.4).  Th('  exitrt'ssions  for  tin'  nonzero  terms  in  tlie  mass  and 
stiffness  matrict-s  given  hy  ('([nations  (8)  throtigh  (15)  can  now  he  exjtrt's.sc'd  in  tc'rms  of  the  matrix  fnnetions 
as  follows; 


where 


=  SfApA) 


Mlj  =  Slj{pA)  +  S}^  {pi  ZZ) 

M;‘j  =  X(i)\{j)[sjj{pA)  +  Sl^iplyy 

Ml,  =  Sfjplxx) 

Klj  =  Sl^iEA) 


K]'.  =  Sj^^iEIzz) 


•i.j 

'tC 


-‘■j 


Kj  =  X{l)X{j)SlP^{EIyy 
Klj  =  50;(G./) 


Xii) 


r  -1  (i  =  2  or  4)  I 
\  1  (Otherwise)  J 


Explicit  ex{)ressions  for  the  five  sets  of  matrix  functions  for  t  =  l.oc  and  j  —  i.oc  are  given  in 
Expressions  for  j  =  1.  /  —  1  are  not  given  because  the  matrix  functions  are  symmetric  (i.e..  5,j  = 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

appendix  A. 

Sj.,)- 


Element  Mass  and  Stiffness  Matrices 

The  nonzero  mass  and  stiffiutss  terms  given  by  equations  (51)  through  (54)  and  equations  (55)  through 
(58).  respectively,  must  be  apjtropriately  a.ssembled  to  form  the  element  mass  and  stiffness  matrices.  This 
procf'dnre  depends  on  the  arrangement  of  the  discrete  degrees  of  freedom  in  the  element  vector  of  unknowns. 
The  first  12  degrees  of  freedom  in  the  vector  of  unknowns  are  the  external  degrees  of  fret'dom  associated  with 
the  cla.ssical  beam  element.  The  higher  order  (internal)  degrees  of  freedom  are  positioned  after  the  external 
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(li'Uii'i’s  of  frccdoiii  in  riic  vector  of  unknovvtis.  S[ieeifi<'ally.  ail  tlie  internal  r/“‘.s  are  gron[>e<|  togetlier.  tlien  all 
the  internal  r/'  s  are  grouped  together,  and  so  ('ii.  The  final  arrangement  of  unknowns  is  givtm  as  follows: 


Mass  Matrix 

Th('  ('lenient  ina.^s  matrix  is  partitioni'd  into  si'veral  suhmat rices  consistent  with  tlu'  arrangc'iiK'nt  of  the 
degrf'cs  of  frei'dom  in  the  vector  of  miknowns  given  in  ('(luation  ((>()).  This  leads  to  a  matrix  of  the  form 

■Mee  m:i  Mb  Mb' 

M}';  0  f)  0 

M=  M'j,  ()  0  (til) 

.‘■(//n  M.'1'i  0 

. 

The  snhniatrix  Mf;^;  is  tt'rmed  the  (’Xternal  snbmalrix  hecaiise  it  is  associatc'd  only  with  the  ('xtc'rnal  degri'es 
of  frf'cdom.  This  siihmatrix.  which  is  given  by  the  following  ecpiation.  is  a  12  by  12  matrix  and  n-ducc's  to  the 
cla.ssical  consistent  mass  matrix  for  constant  crosii-sectional  [iroperties: 
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Tlu'  suhiiMrriccs  Me/'  and  M^y  coujj'.o  the  oxtornal  and  internal  degrees  of  freedom  in  the 

element  mass  matrix.  These  suhmatrices.  whi<h  are  shown  in  ecjnations  (63)  through  (66).  are  variable  sized 
arrays  of  an  order  ecjual  to  I'i  by  the  number  of  internal  degrees  of  freetlom  eorresponding  to  the  partieular 
rontinuous  variable.  For  example,  let  Pn  represent  tlu*  number  of  discrete  degrees  of  freedom  associated  with 
II.  Then,  because  there  are  always  two  external  discrete  degrees  of  freedom  associated  with  n.  the  dimensions 
of  are  1'2  by  P„  —  2.  It  is  of  int('rest  to  note  that  the  coupling  terms  associated  with  the  higher  order 

internal  degrees  of  freedom  become  zero  if  a  sufficient  number  of  internal  degrees  of  freedom  arc'  includc'd. 
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Till'  mass  suhiiiatrict’s  corri’spdiuliiig  to  the  intrriial  dcgm's  of  freedom  are  given  by  the  following  ('(inations 
(e<js,  (f)7)  through  (70)).  Note  that  these  siibinatriees  are  syiniiM'tric  and  banded.  Their  size  depends  on  thi' 
number  of  internal  di'grees  of  freedom  (Miiploved  for  a  particular  eoiMmious  variable.  Given  that  I\  is  thi' 
number  of  disereti'  degree's  of  freedom  associated  with  r  and  that  there  are  always  four  I'xternal  degrees  of 
freedom  associated  with  c.  the  size  of  is  P,  —  -4  by  P,  —  -4. 
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Stiffness  Matrix 


The  stiffness  matrix  is  partitioned  similar  to  the  mass  matrix  and  has  the  form 
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rii('  (’xtcrnal  stifi'iu'ss  suhinatrix  K/r£  is  given  by  the  12  by  12  matrix  in  tlie  following  (‘(jnation  and  redtices 
to  the  classic  >1  bt-am  eletnent  stifftiess  matrix  for  constant  cross-sect ional  ()rop('rt ies. 
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The  submatrices  atid  couple  the  external  and  internal  degrees  of  freedom  in  the  element 

stiffness  matrix  and  are  giveti  by  the  following  equations.  These  submatrice«  have  the  same  dimensions  and 
[troperties  as  the  corresponding  external-internal  coupling  mass  submatrices  (eqs.  (63)  through  (66))  discussed 
Itreviouslv. 
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Tht'  stitfiH'ss  submatrut's  associated  with  the  internal  ilegrees  of  freedom  have  similar  characteristics  as 
theii-  mass  counterparts  and  are  defined  as  follows; 
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All  the  relations  needed  to  completely  define  the  tapered  p-version  beam  element  for  an  arbitrary  number 
of  internal  degrees  of  freedom  are  now  in  hand.  These  include  the  matrix  definitions  given  in  equations  (61) 
through  (80).  the  definition  of  the  cross-sectional  property  polynomials  given  by  equations  (44)  and  (45),  and 
the  ex[)licit  expressions  for  the  five  sets  of  matrix  functions  given  in  appendix  A. 

Numerical  Validation  and  Preliminary  Performance  Analysis 

The  beam  element  developed  herein  is  capable  of  emulating  four  different  types  of  beam  elements:  uniform 
h-version.  uniform  p-version,  tapered  h-version,  and  tapered  p-version.  A  uniform  (h-  or  p-version)  element 
is  created  by  restricting  a\  through  to  be  zero  in  the  generic  cros.s-sectional  property  polynomials.  An 
h-version  (uniform  or  tapered)  element  is  created  by  restricting  Pu  =  Pg  =  2  and  Px<  =  P,„  =  4.  Once  the 
uniform  versions  of  the  element  are  validated,  they  can  be  used  to  approximate  a  tapered  geometry  for  the 
f)urpose  of  validating  the  tapered  versions  of  the  element.  Similarly,  the  p-versions  of  the  element  should 
converge  to  the  same  results  as  the  h- versions  of  the  element. 
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Computational  Approach 

Nuiucrical  ri'sults  Inr  most  of  tht'  casi's  iiiialyzcil  wore  comptitcd  on  ;ui  IBM  Model  80  coiiipiitf’i'  operating 
at  10  MH/  and  e(|nipi)('d  with  an  S03S7  iiuith  coprocessor  and  running  I’t'-DOS  3.3.  Becans('  of  computer 
core  limitations,  those  uniform  h-\’ersion  cases  containing  15  or  more  elements  were  run  on  a  DEC  320(1 
\.\.\  workstation  running  \'MS.  However,  all  computations  used  ('ight-hyte  ffotiting-poiut  r(>pr('sentation  of 
real  numlx'rs  coufortiiing  to  tht'  IEEE,  standard  (tlefault  on  the  IB.M  computer  hut  retiuirt's  a  /Cl  compiler 
ojttion  on  tht'  \’.\\  lomputer).  Bt'catise  tht'  ctunputer  code  used  to  generatt'  tlit'se  rt'sults  was  not  written  to 
take  atlvantagt'  of  tht'  symmetry  anti  hantlt'tlness  of  tht'  global  mass  anil  stiffness  matrici's.  a  fair  cttinparisttn 
of  computt'  tinit's  cannot  bt'  given.  Tht'  only  measure  of  coniputatioiiiil  efficit'iicv  presented  is  the  nutnber  of 
ilegrt't's  t'f  frt't'dom  versus  t  lit'  accuracy  t)f  t  lit'  mtaiels.  Tht'  eigenvalut'  ext  fact  itui  tt'clmitpn'  used  is  a  ctunbiiiat  iou 
of  the  ileterminant  st'arch  method  and  im't'rse  itt'ration  (ref.  Ki).  Tht'  ili'tt'rminant  search  jthast'  i'm;)lt)ys  both 
the  bisection  mt'thotl  iuiil  tht'  secant  mt'thtui  tti  ctuiverge  to  an  eigt'iividue.  Invi'rsi'  iteration  is  tht'ii  ust'ii  to 
obtain  the  associatt'il  eigi'm’i'ctor. 


\  alidation  of  I  niform  h-  and  p-\'ersions 

A  uniiorm  beam  elemt'iit  can  bt'  valitlati'd  rathi'r  t'asily  bt'cause  closeil-form  solutitms  tti  tht'  axial  aiui 
torsituial  vibrations  anti  ‘'nearly  clost'd-fortn"  solutions  to  the  bending  vibratioti  prtibli'in  exist.  (St'e  ref.  12.  ) 
Tht'  ilt'rivatitm  of  the  cltist'd-ftirtn  t'xprossions  fttr  thest'  frt'ijuencii's  is  prt'si'iitt'il  in  appenilix  If  because  tht' 
derivations  comnitinly  fountl  in  tht'  litt'raturi'  lio  not  include  tht-  rotary  int'rtia  ti'rnis. 

The  uniform,  cant ilt'vt'it'd  stt'el  bt'am  with  a  circular  cross  si'ctitui  shown  in  figuri'  3  was  used  tti  valiihitt' 
rht'  uniform  vi'rsions  of  tlu'  t'lt'iiit'nt  ilt'vt'loped  hert'in.  Rt'itresi'iitativt'  h-  aiul  i)-vt'rsioti  finitt'-i'lement  mtiiiels  of 
the  bt'am  art'  ;ilst)  slmwn  in  figurt'  3.  Tht'  matt'rial  [)ropt'rties  of  tht'  bt'am  usi'il  in  the  numt'ricid  stutlii's  wt'ri' 
as  fttllows: 

E  =  30.0  X  10*'  Ib/in'^ 


G  =  11.0  X  10**  Ib/in*^ 


l> 


0.2«d 

380.4 


Ib-sec'^/in'* 


Eor  validation  of  tht'  uniform  vt'rsions  of  the  finite  elt'inent  tievelopi'tl  herein,  the  Itiwt'st  fretpiencit's  of  (his 
beam  wt're  calcttlateii  "t'xactly''  using  the  closed-form  solutions  in  appeiuiix  B  anil  numerically  using  both  tht' 
h-  and  p- version  t'lemt'iits.  Tht'  first  four  bentling  fretpiencies.  tht'  first  torsiotial  fretpiency.  anti  the  first  axial 
frt'tiut'iicy.  from  appenilix  B.  art' 

/,.l  =  27.4051  Hz 


fr2  =  171.0155  Hz 


/l  it  =  478.2()7()  Hz 


frA  =  930.8451  Hz 


fni  =  523.4749  Hz 


/„l  =  841.7997  Hz 

These  six  freijui'iicies  will  bt'  ri'gartleil  as  "exact"  for  the  purjtose  tif  valitlating  thi'  uniform  h-  anti  p- vt'rsions 
of  tht'  beam  eli'int'nt.  Table  I  shows  the  fre(}uencies  obtaineti  using  tlu'  uniform  h-version  t)f  tht'  finitt'  I'li'ini'iit 
as  the  number  of  t'lt'int'uts  tiseil  was  varietl  from  1  to  14.  The  exact  fretiuencit's  are  given  at  tht'  bottom  of 
till'  tablt'  for  rt'ft'renct'.  Although  only  thi'  first  bentling  fretpit'ncy  truly  convt'rgt'tl  with  14  t'li'mi'uts.  tht'  otlu'r 
fretjuencit's  ;irt'  nearly  convt'rgeil  (within  0.00  percent). 

Tablt'  H  shows  thi'  fretiut'ncii's  obtaineti  using  one  p-vt'rsion  elt'inent  as  the  numbt'r  of  intt'rnal  tii'gri'es  of 
frei'iiom  is  increast'd  from  0  to  13  for  t'ach  of  tht*  four  continuous  tlisplacements.  Tluit  is.  in  tablt'  11.  13  inti'inal 
degrees  of  fret'tlom  rt'ft'rs  to  13  internal  tlegrees  of  freeiloin  each  for  i/.  c.  ic,  anti  (L  for  an  aggrt'gtOt'  tot;d  of 
52  intt'rnal  tli'gret's  of  fri'etlom.  Again,  the  exact  fretiuencies  are  noteti  at  tht'  bottom  of  the  tablt'  for  refert'nct'. 
It  is  clear  that  all  thi'  listi'ti  fri'tiui'ncies  have  convergetl  with  13  internal  tli'gret's  tif  freetlom.  anti  that  thi' 
cotivi'rgeil  valut's  are  I'ssi'ntially  t'liual  to  the  exact  values.  In  figures  4  through  9  thi'  information  in  ttdili's  1 
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and  II  is  plotted  in  terms  of  the  percent  error  versus  tiie  number  of  degrees  of  freedom.  The  plots  clearly  show 
that  both  tilt'  p-\ersioii  element  and  the  h-version  element  converge  to  the  exact  answer,  and  that  the  p- version 
elt'nu'iit  has  a  faster  convergence  rate  than  tlit'  h-version  element. 

Validation  of  Tapered  h-  and  p- Versions 

Tht'  tapt'rt'd  steel  beam  with  circular  cro.ss  section  shown  in  figure  10  was  used  to  validatt'  the  tapert'd 
\'ersions  of  the  finite  ('lenient  developed  herein.  Also  shown  in  figure  10  are  represt'iitative  h-  and  p-vt'rsion 
models.  Because  the  rt'sults  shown  in  tables  I  and  II  have  validated  the  uniform  versions  of  the  beam  eh'inent. 
till'  exact  fretiuencies  of  the  tapered  tieani  can  be  determined  with  a  uniform  h-version  model  if  a  sufficiently 
large  number  of  eleiiu'iits  are  employt'd. 

Tht'  fretiuencit's  obtained  for  the  tapt'rt'd  beam  using  uniform  h-version  finitt'  t'lenients  art'  prt'st'iitt'd  in 
tablt'  III.  It  should  bt'  noted  that  tht'  fretjut'iicies  in  table  III  approach  the  converged  values  from  below 
(whereas  in  tables  I  and  II.  the  frt'qnencies  converge  from  above).  Also,  the  number  of  elements  retinired  to 
rt'ach  convt'rgence  for  the  tapered  beam  is  over  seven  times  that  rerjuirt'd  for  a  comparablt'  uniform  bt'ani. 
Both  tht'st'  behaviors  can  be  explained  by  the  fact  that,  as  uniform  elements  are  added,  the  geometry  of  the 
inodt'l  is  changing.  Tlie  geomt'try  approaches  a  smoothly  tapered  beam  as  the  number  of  elements  in  the  model 
approacht's  inhnity.  Here,  the  model  is  assumt'd  to  bt'  converged  with  100  elt'inents  (000  degrt't's  of  frt'edom). 

Tht'  frt'tint'ticit's  prt'dicted  by  using  a  tapt'red  h-version  finite-element  model  are  shown  in  table  I\'.  In  this 
cast'  the  gt'omt'try  of  the  problem  is  represented  exactly  using  one  element,  and  therefore  the  geometry  of  tht' 
modt'l  is  not  changed  as  the  number  of  elements  is  increast'd.  It  should  be  noted  that  now  convergence  is  from 
above  and  is  achieved  using  approximately  one  seventh  as  many  elements  as  in  table  III. 

Table  \'  shows  the  frecjnencies  which  were  obtained  fttr  the  tapered  beam  using  one  tapered  p-version 
element  and  varying  tlie  number  of  internal  degrees  of  freedom  from  0  to  19.  With  19  internal  degrees  of 
freedom.  22nd  order  polynomials  are  used  for  lateral  bending  and  20th  order  polynomials  for  the  axial  and 
torsional  displacements.  The  rea.son  for  showing  the  results  of  .such  a  high  order  element  is  to  demonstrate  the 
numerical  stability  of  the  element.  It  is  seen  that  the  p-version  results  in  table  V  converge  very  close  to  tlie 
values  prediett'd  by  the  h-version  beams  in  tables  III  and  IV  but  with  the  use  of  far  fewer  degrees  of  freedom. 

The  data  in  tables  III  through  V  are  graphically  depicted  in  figures  11  through  10  in  terms  of  percent  error 
versus  the  number  of  dt'grees  of  freetlom.  These  figures  show  the  dramatic  improvement  in  convergence  which 
is  rt'alized  by  using  tapered  h-  and  p-version  elements  in.stead  of  uniform  h-version  elements. 

Conclusions 

The  derivation,  and  validation,  of  a  new.  tapered,  p-version  beam  element  which  both  facilitates  convergence 
checks  and  produces  a  better  convergence  rate  than  nontapered.  h-version  beam  elements  has  been  described. 
Tht'se  two  characteristics  complement  each  other  and,  when  combined,  provide  a  powerful  and  versatile  beam 
('lenient  which  is  easy  to  use.  The  shape  functions  on  which  the  element  is  based  were  derived  by  using 
orthogonality  rt'lationships  w'hich  produce  element  matrices  that  are  extremely  well-conditioned  and  of  a  form 
allowing  explicit  integration  in  the  derivation  of  the  element  matrices.  The  latter  feature  eliminatt's  the  need 
for  numerical  quadrature:  thus,  roundoff  error  is  reduced.  The  shape  functions  are  hierarchical  such  that  higher 
order  element  matrices  can  use  the  element  matrices  from  previous  lower  order  analyses.  This  simplifies  the 
derivation,  coding,  and  validation  of  the  element.  The  present  form  of  the  beam  element  has  been  derived  in  a 
manner  which  allows  for  an  infinite  number  of  intt'rnal  degrees  of  freedom.  The  bt'am  element  has  bt'en  tt'stt'd 
with  up  to  22nd  ordt'r  C^-type  shape  functions  and  up  to  20th  order  C^’-type  shape  functions  with  no  evidence 
of  ill  conditioning  or  significant  roundoff  error. 

N.^SA  Langley  Resparcli  Center 
Hain[)t()n.  VA  2.'566o-.j22.'j 
.June  23.  1989 
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Appendix  A 

Expressions  for  Matrix  Functions 

Explicit  ('xpix'ssions  for  the  matrix  Functions  are  preMuiteci  in  tliis  appendix.  Each  matrix  function  is  a 
function  of  the  ehmient  len‘>;tli  /  and  the  coidficituits  from  tfie  geiu'ric  cross-sectional  j)ro[)erty  (jolvnomials.  u,) 
through  a  (. 

The  df'finition  of  the  5*'  set  is 


(Al) 

(A  2) 

(A.T 

(A4) 


5'i’,  =  -  3  [(/  -  -  685/'^  +  4990/-  -  14910/  +  15492)//.}/-^ 

+  2(/  -  8)(2:j/"‘  -  41G/'^  +  2097/^  -  7419/  + 

+  (j(i  -  8)(/  -  7)(8/'^  -  90/2  ^  322,  _  3G5)fi2/^ 

+  4(/  -  8)(/  -  7)(/  -  6)(/  -  4)(4/  -  5)fJi/ 

-h  'Aoli  -  8)(/  -  7)(/  -  6)(/  -  5)(2/  -  r))ao] 

/ [1(38(587/""  -  15825/'^  +  1661  lo/'^  -  846075/2  -p  2087798/  -  1994340)]  (A5) 


5'/,,  =0  (A6) 

=9. St 

S^l,  =  /v/2/  -  3  [{283/""  -  7825/'*  +  84335/^  -  442025/2  -f-  1124532/  -  1 108260 )f/4/^ 

-f  5(/  -  8){50/"*  -  995/^  +  7224/2  _  22641/  -f  25770)a:i/-^ 

+  3(/  -  8)(/  -  7)(68/^  -  905/2  -f-  3907/  -  544())«.2/2 
+  3(/  -  8){/  -  7)(/  -  6)(52/2  -  413/  +  785)«i/ 

+  35{/  -  8)(/  -  7)(/  -  6)(/  -  5)(4/  -  13)ao] 

/[168(587/"  -  15825/"*  +  166115/"’  -  846075/2  -p  2087798/'  -  1994340)]  (A7) 

5.^',  =0  (A8) 


5,",  =  l{2i  -  3)  [3(7/'*  -  42/'*  -+•  7/2  +  168/  +  60)a4/'* 

+  (2/  +  3)(2/  -  9){7/'2  -  21/  -  10)«;!f'‘ 

-P  2(2/  -P  3)(2/  -  9)(5/:2  -  15/  - 
+  4(2/  +  3)(2/;  +  1)(2/:  -  7)(2/  -  9)(rti/  +  2«())] 

/[16(2/  +  3)(2/:  +  1)(2/:  -  1)(2/;  -  3)(2/  -  5){2/  -  7)(2/  -  9)]  ( A9) 
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=  -/- v'27^v/27^[{i(/  -  l)(r  -  2/  -  5)fi4/‘* 

-f-  {.'  -  1  }(7r  -  14/  -  :ib}uj/’  +  2(/  -  1)(2/  +  3)(2/  -  7)(«2^  +  "l )] 

/  [U)i2/  +  3)(2/  +  1 1(2/  -  1)(2/  -  3)(2;  -  5)(2/  -  7)]  (AlO) 

=  -/v'2/  -  :W2i  +  l[l3(r  -  /  ~  1())(/'^  -  /  -  8)a4/^ 

+  s(/  +  2)(/  -  3)(2/  +  3){2/  -  +  H(2i  +  5)(2/  -  7)(2/'^  -  2/  -  9)a2*^ 

+  8(2/  +  5)(2/  +  3)(2/  -  5)(2/  -  7)(«i/  +  2(i())j 

/[fiU2/  +  5)(2/  +  31(2/  +  1)(2/  -  1)(2/  -  3){2/  -  5)(2/  -  7)]  (All) 

=  /\'2/  -  3v/2/  +  3[l()/(/'-  -  7)/i4/'*  +  /(13/''-’  -  8r))«:j/- 
4-4/(2/  +  ■>)(2/  -  oKrto/  +  «i )] 

/[32(2/  +  5)(2/  +  3)(2/  +  1)(2/  -  1)(2/  -  3)(2/  -  5)]  (.412) 

=  -/V2/ -3\/2/  +  3[(/'  +  1)/(13/''^  +  13/  -  llG)/i4/'-^ 

+  (  /  +  l)/(2/  +  7)(2/  -  5)(3(i:i/  +  2«2)] 

/[32(2/  +  7)(2/  +5)(2/  +  3)(2/  +  1)(2/  -  1)(2/  -  3)(2/  -  5)]  (A13) 


,()  / * \/2 /  —  3 \/2 i  +  I  (i  +  2)(/  +  l)i{2u.il  + //,j ) 

"  32(2/  +  7)(2/  +  5)(2/  +  3)(2/-  +  1)(2/  -  1)(2/  -  3) 


(AM) 


^.,(1  —  \/2i  —  3v^2/  +  9(/  +  3)( /  4-  2){/  4- 

"  64(2/  +  9)(2/  +  7)(2/'  4-  5)(2/  +  3)(2/  +  1)(2/-  -  1)(2/  -  3) 


(A15) 


The  definition  of  the  5*'^  set  is 

and  till’  e.xplicit  integrations  ar/' 

12«4/‘  4-  lofi;,/^  4-  2f)c2l'^  +  3()rti/  4-  6()no 

‘Vl  ”  60/ 

(1^7  T2/i4/^  4"  lofi.j/'^  +  20(121^  4-  3()//i/  4-  60//() 

LJ  ""  60/ 

r2//4/’  4-  laa.j/'^  4-  20a2l^  +  OOrt]/  4-  60//() 

2.2  ""  60/ 

=  V2i  -  3[3(2<J/''’  -  406/:’^  +  1833/  -  2644)//4/'^ 

'  =  4-  3(/  -  6)(25/^  -  202/  4-  387)//;!/'^  4-  2(i  -  6)(/  -  r))(26/  -  83)//2/ 

+  10(/  -  6)(/  -  r>)(/'  -  4)/i|]/[30(47/'*  -  628/'^  +  2729/  -  3792)] 


(AlO) 


(A17) 


(A  18) 


(A19) 


{A20) 


(A21) 
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.s-'l'.'’  =0  (A22) 


=  -V2)  -  ;5[3r2yr*  -  wtir  +  lb:5:5/  -  2(U4)a.jr*  +  :H/  -  (j)(25;^  -  202/  +  387)«:j/- 
,  =  ^  2(1  _  (j){;  _  5)(2(i/  —  83)m2/  +  !(•('  “  0)(/  —  5)(/  —  4)rti] 

/  [;50( 47/'*  -  028/-  +  2729/  -  3792 )]  ( A23 ) 


S]'l'  =  'pbi^  -  210/'*  +  305/-  +  30/  -  88)(2/  -  3)//.,/' 

'  =  '*-^  +  2(2/  +  1){2/  -  3){2/  -  7)(5/'-^  ~  15/  +  TUr/' 

+  4(2/  +  l)(2/  -  31(2/  -  7)(3/'^  -  9/  -  4)/)2/- 

-r  4(2/  +  1  ){2/  —  1 1(2/  —  3)( 2/  —  •))(2/  —  7 )(// j /  +  2//|) )j 

/[S/(2/  +  11(2/  -  11(2/  -  31(2/  -  51(2/  -  7)] 

Vi'fli  =  -y27^/27^[2(/  -  1)(7/'^  -  14/  -  8)//.,/'^ 

+  3(/  -  1)(5/-  -  10/  -  ())//;j/'^ 

-t-  4 ( /  —  1 ) ( "2/  +  11(2/  ~  5)( //'/^  4"  1 )] 

/[8(2/  +  11(2/  -  11(2/  -  31(2/  -  5)] 

.sT.flj  =  v^2/'  -  :W2>  +  1(21(/  -  r)(7/''^  -  7/  -  2G)//4/'* 

+  /(/■  -  1 1(2/  +  3){2/  -  5)(3//;(/'^  +  2//2/)] 

/[8(2/  +  3)(2/  +  11(2/  -  11(2/  -  3)(2/  -  5)] 


y/^  3v/27T3(/  +  l)/(/  -  1)(2//4/^  +  r/3/‘) 
8(2/  +  3 )  { 2  /  +  1 )  ( 2  /  —  1 )  ( 2  /  —  3 1 

\/2/  —  3\/'2i  +  5{/  +  2)(/  +  l)/(/  —  1  j/ij/'* 

10(2/  +  5 )  ( 2  /  +  3 )  ( 2  /  +  1 )  ( 2  /  —  1 )  ( 2 1  “■  3 ) 


oi)P 

•-  i.i+j 

/  =  4.x.  j=.j.x 


=  0 


'rh(‘  (h'finitioii  of  thr  5*  sot  is 
rl 


Slj  =  ^  P2(.O.V‘A'j  dr 


(A24) 


(A25) 


(A20) 


(A27) 

(A28) 


(A '291 


(A30) 


(A31) 


I 
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aiui  tlu'  oxplk'it  intei^rations  art' 


/( +  ->4(^1/  -i-  234(i()) 
.S 1  1  =  - 


4.1 


()30 


/“  1 1  /  (■/■)/“■  -l-  42ti  ]  /  4-  1 42f/() ) 
^  2520 


I  /(liiaj/"  -h  81tJi  /  +  lti2fi()) 

“  r2(j() 


.1  /■(  I'Jf/o/'  +  3(i«i/  -t-  TSfit)) 


■'’i.t  -  ~ 


2520 


I  I'H  -iii  >l~  +  !)« 1  /  +  24f;()  ) 
=  ^4520 


/'(25r)_)/“  -I-  12a  [/  -h  7^a(|) 
""  2520 


^,|  /'^( at/'j/"  +  !)a  1  /  -f-  ISuo) 

""  ^720 


/( 145a2/‘  +  I'^^Oai/  +  234ri()) 
=  - — - 


030 


,]  +  OOcii/  +  I32a()) 

2520 


4...  -  - 


cl  _ 

•^4.4  - 


/■*(10a2/‘  +  15a,/  +  24a„) 
2520 


(A32) 


(A33) 


1A34) 


:A35) 


A3()) 


A  37, 


( A3,'<) 


(A3!0 


;  A  lO  i 


(All) 


•S’J ,  =  -  77250*  +  !  140200*  -  827175r  +  2047034/  -  4122OOO)a20* 

'  =  "’■*"  +  (/  -  10)(  173/*  -  5005/-*  +  53215/'^  -  245015/  +  410202)//,/^ 

+  2(/  -  10)(/  -  0)(/  -  7)(00/'^  -  1130/  t  3302)ao/] 

/[1200(4750/  ’  -  174100/*  +  2512445/**  -  17844110/^  +  02404800/  -  85080000)]  ( A 42) 

5*,  =0  (A43) 

,  =  1 1 


.S']  ,  =  y/2i  -  5  [{02/'*  -  2315/  *  +  34000/**  -  245305/^  +  808878/  -  r208400)a2/* 

'  =  ■>.10  ^  .2{,  -  10)(/  -  8)(/  -  0)(33/'^  -  408/  +  1400)rti/** 

+  (/  -  10){/ -  0)(20/**  -  001/'^  +  4001/ -  11470)a,)/‘*] 

/[1008()( 434/'’  -  15875/*  +  220000/**  -  1028305/'“*  +  5701080/  -  7801020)]  (All) 

.S'.],  =0  (A  45) 
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Appendix  B 

Derivation  of  Vibrations  of  Uniform  Beam 

The  (l('ri\ati()u  of  tlif  clo.st'ci-foMii  solutions  to  tli<‘  vibrations  of  a  iiniforni  beam  is  pr(>s('ntc(l  in  this  ajtjx'udix. 
I  ii.'i  il.v  bfiidin”  vit)ration  probhun  is  solved,  and  tluni  ttie  torsional  aiui  axial  vibiailon  pix/bl. 

Consider  the  translational  displactniient  r  only.  Csing  the  ('xjin'ssions  for  the  kiiu'tic  and  strain  ('iiergii's 
yivi'ii  in  ('(luations  (21  ;ind  (d).  respect ivi'ly.  and  assnniiu{>;  hannonic  motion  give  the  following  governing 
ditfennitial  equation: 

r""+^c"-^4^e  =  ()  (Bl) 

When  taking  into  account  the  cantiii'ver  conditions,  tlu'  geometric  boundary  conditions  ar<' 

(•(())  =  ()  (B2) 


and  the  mitural  bonndarv  conditions  are 


r'(0)  =  0 


r"(7)  =  0 


t 


The  general  solution  to  ecpiation  (Bl)  is 


r(./')  —  -f-  A*2  c'  -I-  A’;j  cos(/y2.r)  -h  Aq  sin(/>2.r)  (Bli) 

where  Aq.  Aqi.  A'q.  and  Aq  are  cotistants  that  depend  upon  the  bomi<iary  conditions,  and  /q  and  1)2  art'  didiiK'd 
by 


^  A  ^  A-  .4  -2 


EIzz 


(BIO) 


■Substituting  th('  general  solution  given  by  eiination  (B(i)  into  the  boundary  conditions  (eejs.  (B2)  through  (BT))) 
results  in  four  linear  equations  for  Aq  through  A',)  tliat  can  be  expresse<i  in  matrix  form  as  follows; 


hj  l)\  (‘  -62COs(/)2/)  -I)2^\n{h2l)  j  A'.-j  (  ^  0 

.h](hf  -i-  .dq)  c^'i^  +  -U)  c“^''^  ~  d2)sin(/i2/)  ~l>2il>2  ~  -'U ) ''<>s(/r2/) .  I  Aq  J  f  0 


(BID 


The  bending  natnr  d  freciuencies  of  the  beam  are  those  which  make  the  determinant  of  the  1  by  1  matrix 
ill  (’(luaMoii  (Bll)  e(|nal  to  zero.  For  the  example  dc'scrilied  in  the  main  text  of  tins  report,  the  first  four 
frequencies  art' 
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f,l  ^  27.4651  Hz 
/,,2  =  171.6155  Hz 
/,,3  =  478.2670  Hz 
=  930.8451  Hz 

Xcnv  consider  only  the  torsional  displacement  0.  Taking  equations  (2)  and  (3)  as  the  kinetic  and  strain 
energy's,  respectively,  and  assuming  harmonic  motion  giv-e  the  following  governing  differential  equation; 


The  boundary  conditions  are 


0"  + 


plxx^l 

GJ 


6  =  0 


(B12) 


61(0)  =  0 
6'(l)  =  0 

With  h~  =  ^ .  the  general  solution  to  equation  (B12)  can  be  written  as 

6{x)  =  ki  CQs(bx)  +  k2iihi{bx) 


(B13) 

(B14) 


(B15) 


Substituting  equation  (Bio)  into  the  boundary  conditions  leads  to  the  2  by  2  matrix  system  of  linear  equations: 


1  0 
L-sin(6/)  cos(6/) 

Setting  the  determinant  equal  to  zero  and  solving  for  ujq  yields 


/  1  _  /  0  1 
U2  J  "  I  0 / 


ntr  /  GJ 


=  TH-i 


2/  V  P^XX 

For  the  beam  of  figure  3.  the  first  torsional  frequency  is 


(n  =  1.2....) 


(B16) 


(B17) 


fgi  =  523.4749  Hz 

riic  gcjverning  diff<’rential  equation  for  the  axial  displacement  of  the  beaiii  is  of  the  same  form  as  that  for 
the  torsional  displaonient  given  in  equation  (B12).  Only  the  constant  b  is  different.  The  natural  frequencies 
can  be  shown  to  be 


For  the  beam  of  hgnre  3.  the  first  axial  frequency  is 


(n=  1.2....) 


(B18) 


/„l  =  841.7997  Hz 
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ThMo  I  Convergence  of  Frequencies  for  Uniform.  li-V^ersion  Beam 


\(iof 

.\elem 

/,!.  Hz 

fr2-  Hz 

fri-  Hz 

/,,4.  Hz 

fei-  Hz 

ful-  Hz 

6 

1 

27.5953 

270.3419 

577.2137 

928.2170 

12 

2 

27.4784 

173.0589 

581.7178 

1649.2712 

537.0115 

863.5679 

18 

3 

27.4079 

172.1745 

484.1519 

1081.3530 

529.4733 

851.4457 

2-1 

4 

27.4600 

171.8138 

481.9206 

944.1763 

526.8447 

847.2186 

30 

5 

27.4055 

171.7006 

479.9624 

941.5603 

525.6302 

845.2656 

,31) 

ti 

27.4053 

171.0575 

479.1306 

936.7340 

524.9711 

844.2057 

12 

7 

27.4052 

171.0385 

478.7476 

934.2196 

524.5739 

843.5670 

■18 

8 

27.4052 

171.0291 

478.5541 

932.8928 

524.3162 

843.1526 

54 

y 

27.4652 

171.6240 

478.4485 

932.1524 

524.1396 

842.8685 

(10 

10 

27.4052 

171.0211 

478.3871 

931.7163 

524.0132 

842.6654 

titi 

11 

27.4652 

171.6193 

478.3496 

931.4469 

523.9198 

842.5151 

.  72 

12 

27.4052 

171.6182 

478.3256 

931.2736 

523.8487 

842.4008 

78 

13 

27.4052 

171.0175 

478.3097 

931.1583 

523.7934 

842.3119 

84 

14 

27.4652 

171.6170 

478.2988 

931.0791 

523.7495 

842.2413 

Exact  .  .  . 

27.4051 

171.6155 

478.2670 

930.8451 

523.4749 

Tai)le  II.  Convergence  of  Frequencies  for  Uniform,  p- Version  Beam 


Ndof 

Idof 

fri.  Hz 

/c2-  Hz 

fvi-  Hz 

/c4.  Hz 

6 

0 

27.5953 

270.3419 

577.2137 

928.2170 

10 

1 

27.4734 

173.1390 

906.1908 

844.9599 

14 

2 

27.4052 

r2.5710 

490.8033 

841.8572 

18 

3 

27.4651 

171.0203 

490.0998 

22 

4 

27.4651 

171.6178 

478.4139 

987.4482 

523.4749 

841.7997 

20 

5 

27.4651 

171.6155 

478.4039 

841.7997 

30 

0 

27.4651 

171.6155 

478.2675 

841.7997 

34 

7 

27.4651 

171.6155 

478.2674 

841.7997 

38 

8 

27.4651 

171.6155 

478.2669 

523.4749 

841.7997 

42 

9 

27.4651 

171.6155 

478.2669 

523.4749 

841.7997 

10 

10 

27.4051 

171.6155 

478.2669 

523.4749 

841.7997 

50 

11 

27.4651 

171.6155 

478.2669 

930.8451 

523.4749 

841.7997 

54 

12 

27.4651 

171.6155 

478.2669 

523.4749 

841.7997 

58 

13 

27.4651 

171.6155 

478.2669 

523.4749 

841.7997 

Exact  .  .  . 

27.4651 

171.6155 

478.2670 

930.8451 

523.4749 

841.7997 
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lahlc  III.  C’i'iiV('rg('n(  t'  of  Frcciucm  irs  for  Modt'liiii!,  Tapered  Beam  With  riiifonii  Beam  Klemeiits 


'  N(h)f 

1 

Xelem 

/,  (.  Hx 

/,  >.  Hz 

Hz 

/r4.  Hz 

hi-  Hz 

fnl-  Hz 

ti 

1 

30.7113 

358.0025 

577.2137 

928.2170 

12 

2 

30.t)302 

172.. 3000 

709.7730 

1391.3137 

934.4910 

1255.9180 

IS 

3 

71.0803 

187.9,347 

455.8570 

1109.0279 

1089.9504 

1,5J8.1870 

2  4 

4 

77.8001 

205.9470 

441.9430 

893.5114 

1150.4921 

1381., 3930 

3(1 

5 

80.0403 

218.2811 

455.8057 

834.3421 

1189.1758 

1390.3582 

3G 

0 

82.0833 

220.0285 

470.7080 

8,35.1007 

1207,3211 

1404., 304  5 

;  42 

7 

83.7385 

230.9900 

481.8723 

848.0075 

1218. .3055 

1409. 1.399  . 

!  4S 

84.4037 

234.3112 

489.7455 

800.5784 

1225.5007 

1412,2157  ! 

34 

0 

84.0334 

230,0301 

495.3418 

870.3900 

1230.5101 

1414.3128 

bO 

10 

83.3083 

238.3102 

499.4259 

877.7242 

1234.0013 

1415.8008 

()G 

11 

83.3707 

230. .3054 

502.4900 

883.2409 

1230.0807 

1410,9088 

!  72 

12 

fill 

240.5273 

504.8484 

887.4925 

1238.0840 

1417,74,50  ! 

7S 

13 

83.0270 

241.2805 

500.7015 

890.8284 

1240.2401 

1  118.3940 

14 

80.0520 

241.8811 

508.1840 

893.501 1 

1241.4740 

1418.9092  ! 

'JO 

13 

80. 1520 

2  42.3070 

509.3901 

895.0784 

1242.4707 

1419.3238 

;  130 

23 

80.5808 

244.4943 

514.7104 

905.4243 

1240.7811 

1421.1123 

:  300 

30 

80.7753 

245.4031 

517.0280 

909.7027 

1248.0001 

1421.8042 

430 

7") 

80.8097 

245.5723 

517.4313 

910. .5840 

1248.9370 

1422.0033 

!  000 

100 

80.8373 

245.7077 

517.8087 

911.2405 

1249.2228 

1422.1210 

Table  I\’.  Couvcn'^enct'  of  Fretjueiieies  for  Tapert'd.  h-\’ersioii  Bc'iim 


Xdof 

.X'eiem 

fr\-  Hz 

/rz.  Hz 

fr:i-  Hz 

/, ).  Hz 

(i 

1 

80.9704 

207.8185 

1273.8317 

1434.05,50 

12 

2 

80.9005 

251.7443 

,533.1733 

1203.8027 

1271.8282 

1432.0379 

18 

3 

80.8,501 

247.4,589 

538.74.30 

940.3004 

1203.2850 

1428.3205 

24 

4 

80.8434 

240.2835 

520.0757 

958.9185 

12.58.0545 

1425.9792 

30 

•*) 

80.8390 

245.9.340 

,521.0733 

934.7022 

1255.1034 

1424.7097 

.30 

0 

80.8384 

24,5.8121 

519.2780 

921.7052 

1253.4.584 

1423.9058 

42 

7 

80.8379 

245.7028 

518. .5013 

910.4307 

12,52.3820 

1423.4975 

48 

8 

80.8,370 

245.7400 

518.2391 

914.1021 

1251.0047 

1423.1850 

54 

9 

80.8375 

245.728.3 

518.0774 

912.9740 

1251.1033 

1422.9009 

00 

10 

80.8374 

245.7210 

517.9880 

912.3740 

1250.8000 

1422.8087 

00 

11 

80.8,374 

245.7170 

.517.9345 

912.0208 

12,50.5280 

1422.()900 

72 

12 

80.8,374 

245.7149 

517.9004 

911.8101 

12,50.3208 

1422.(i002 

78 

1.3 

80.8373 

245.7131 

.517.8775 

911.0072 

12,50.1582 

1422.5294 

84 

1 4 

80.8,373 

245.7119 

.517.8015 

911., 5087 

1250.0287 

1422.4729 

Tat)l(’  Convorgenco  of  Froquoncios  for  Tajx'rcd.  p-V^ersioii  Beam 


Ndof 

Idof 

frl.  Hz 

fv2-  Hz 

/,.3-  Hz 

fr4.  Hz 

foi-  Hz 

ful.  Hz 

0 

80.9704 

207.8185 

1273.8317 

1434.6550 

10 

1 

80.9555 

249.0871 

007.0043 

1254.6750 

1434.2370 

14 

2 

80.8075 

248.0820 

545.9870 

1173.2573 

1249.4274 

1422.2183 

18 

80.8417 

240.2815 

537.1999 

1018.3391 

1249.2122 

1422.1270 

22 

4 

80.8376 

245.7577 

520.7508 

981.2253 

1249.2122 

1422.1223 

20 

80.8374 

245.7151 

517.9835 

920.7022 

1249.2093 

1422.1164 

;i() 

0 

80.8374 

245.7151 

517.9218 

912.1356 

1249.2074 

1422.1149 

34 

7 

80.8373 

245.7127 

517.8910 

912.1356 

1249.2067 

1422.1145 

38 

8 

80.8373 

245.7101 

517.8475 

911.6613 

1249.2065 

1422.1145 

42 

9 

80.8372 

245.7087 

517.8238 

911.3808 

1249.2064 

1422.1144 

40 

10 

80.8.372 

245.7080 

517.8140 

911.2864 

1249.2064 

1422.1144 

50 

11 

80.8372 

245.7078 

517.8104 

911.2582 

1249.2064 

1422.1144 

•"j4 

12 

80.8372 

245.7077 

517.8092 

911.2498 

1249.2064 

1422.1144 

58 

13 

80.8372 

245.7077 

517.8088 

911.2473 

1249.2064 

1422.1144 

02 

14 

80.8372 

245.7077 

517.8087 

911.2466 

1249.2064 

1422.1144 

00 

15 

80.8372 

245.7077 

517.8086 

911.2464 

1249.2064 

1422.1144 

70 

10 

80.8372 

245.7077 

517.8086 

911.2463 

1249.2064 

1422.1144 

74 

17 

80.8372 

245.7077 

517.8086 

911.2463 

1249.2064 

1422.1144 

78 

18 

80.8372 

245.7077 

517.8086 

911.2463 

1249.2064 

1422.1144 

82 

19 

80.8372 

245.7077 

517.8086 

911.2463 

1249.2064 

1422.1144 

34 


00 


p-version 

h-version 


”ur<‘  •) 


i^urc 


Number  of  degrees  of  freedom 

Pcrci'iu  ('’Tor  v('rsus  uumlx'r  of  (li'f^n'cs  of  fnM'doin  for  socoiui  hi'iidins  frcHjtu'ncv  of  uniforiii  licam 


Percent  error 


Number  of  degrees  of  freedom 


j.  Percent  error  versus  nuinher  of  degrees  of  freedom  for  t  hird  hending  frecinency  of  uniform  tx'am. 
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Percent  error 


Number  of  degrees  of  freedom 


Percent  error  vin'Mis  numl)(>r  of  degrees  of  freedom  for  fourtli  Ix'iuiing  fre  (lumicy  of  uniform  Ix'am 


Percent  error 


Number  of  degrees  of  freedom 


Figure  r^rueent  error  versus  numher  of  degrees  of  freedom  for  first  torsional  frrHiuency  of  uniform  beam 


100.0 


Tapered  p-version 
Tapered  h-version 
Uniform  h-version 


P  inure  11.  Percent  I'rror  ver.‘'U.s  luuuher  of  degrees  of  freedom  for  first  bending  frequency  of  tapered  t)eain. 


Percent  error 


figure  12.  Percruit  error  vitsus  nuinlx'r  of  degrees  of  freedmn  for  second  bending  fre<iuency  of  tapi'red  Ixuiin 


40 


Percent  error 


Tapered  p-version 
Taoered  h-version 
Uniform  h-version 


Number  of  degrees  of  freedom 

Fi‘j;ur('  13.  Pt'rci'iit  ('rror  versus  ninnher  of  duffrees  of  freedom  for  third  Ixuiding  fr('(iueiu'y  of  tapc'red  beam. 


Percent  error  .|  q 


—  Tapered  p-version 

—  Tapered  h-version 

—  Uniform  h-version 


Number  of  degrees  of  freedom 


Figure  14.  Percent  errf)r  versus  nunil)er  of  degrees  of  freedom  for  fourtli  t)ending  freciueney  of  tapi'rc'd  Ix'am 


Tapered  p-version 
Tapered  h-version 
Uniform  h-vcrsion 


I-iSurt'  15.  Percent  error  versus  number  of  degrees  of  freedom  for  first  torsional  frequency  of  tapered  beam. 


Figure  16.  Percent  error  versus  number  of  degrees  of  freedom  for  first  axial  frequency  of  tapered  beam. 
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